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SUMMARY 

This reporc la volumtt X of a three volume aeries entitled 
"Free Wake Techniquen for Rotor Aorodynamtc Analysis'* and covering the 
following topics: 

Volume I (present volume) "Summary of Rasults and Background Theory" 
reviews the results obtained to date using both complete and 
simplified wake models and summarizes the theoretical background on 
which these models are based. 

Volume II "V^ortcx Sheet Models" (Ref. 5) presents the results of 
computations using complete and modified vortex sheet models and 
tests the sensitivity of the solutions to various assumptions used 
in the development of the models. The complete codings are included. 

Volume III "Vortex Filament Models" (Ref. 6) discusses results 
obtained using a vortex filament model, as opposed to sheets, again 
using various modelling techniques and including the computer codings. 
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SECTION I INTRODUCTION 

This report will discuss results obtained to date during the developoent 
of a consistent aerodynaalc theory for rotors In hovering flight. Methods of 
aerodynamic analysis have been developed which arc adequate for general design 
purposes until such time as core elaborate solutions become available, in 
particular solutions which include real fluids effects. It Is recognized that 
such solutions arc not yet within the capability of even the most powerful 
computational facilities. Nevertheless it is hoped that continued research, 
both experimental and computational, will lead to further definition of the 
real fluids problem and narrow the range of empiricism now required to handle 
these effects. 

Several problems were encountered in the course of this development, 
and many remain to be solved, however it is felt that a better understanding 
of the aerodynamic phenomena involved has been obtained. Remaining 
uncertainties are discussed in the following sections and in the companion 
volumes (Refs. 5 and 6) covering this study. 

Experimental investigations (Refs. 1, 2, 3, 4) have shown that the 
wake geometry of a hovering rotor differs appreciably from the simple spiral 
form of constant diameter assumed in classical theory. The airload distribution 
on a rotor blade is critically dependent on this geometry and particularly 
on the distance between the blade and the tip vortex generated by the previous 
blade at their first encounter, as indicated in Fig. 1, and on the structure 
of the vortex at this encounter. There is strong experimental evidence that 
this vortex, trailed initially from the blade as the bound circulation drops 
rapidly from a peak at about 90Z of rotor span to the tip, rolls up rapidly 
into a tight and strong tip vortex whose strength is close to the maximum 
strength of the bound circulation on the rotor blade. What is not known is 
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the azimuth locationa over which this roll up occura, the structure of the 
resultant vortex core and the mutually induced velocities in the near wake 
during the roll up process. These velocities are important in determining 
the vertical location of the vortex at first encounter. 

It is convenient to consider the wake as divided into three segments. 

The first segment, the near wake, contains that portion of the wake attached 
to a particular blade and before first encounter with a following blade. The 
second segment, the intermediate wake, includes the wake from the first 
encounter through two or more spirals. The third segment, the far wake, 
contains the rest of the wake to infinity. It is generally a reasonable 
approximation to treat only the first two segments of the wake as free, that 
is with a geometry determined by the induced velocities in the wake, and to 
consider the far wake as a vortex cylinder of constant radius and strength 
determined by the radius and vortex spacing at the end of the Intermediate 
wake. Calculation of the intermediate and far wake geometries requires 
extensive computational capabilities, but otherwise present no serious problem. 
There is some question as to the suitability of using a fixed geometry far 
wake, as discussed in section IV, however, as might be expected, assumptions 
relating to the far wake do not appear to have a major impact on the solutions. 

Of more concern is the treatment of the near wake. Figure 2 shows a 
typical spanwise distribution of bound circulation on a two-bladed rotor from 
both the experimental results of Ref. 4 and the analytical results of Ref. 7. 
The rapid drop off of circulation over the outer portion of the blade results 
in a trailing wake system which may be modelled mathematically as a scries of 
curved vortex filaments as shown in Fig. 3. The induced velocities acting 
on this system of vortex filaments from the far and intermediate wakes can 
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bd readily dccenslned given the wake geometry. However it is primarily the 
mutually induced velocities due to the near wake which determines the roll 
up process and hence the eventual near wake displacements. 

In section It of this report various models are discussed for determining 
the near wake displacements and it is shown that reasonable agreement with 
experimental results arc obtained with versions of the following three models 

a) A model in which the near and incevmedlate wakes are 
represented by spiral vortex sheets which are assumed 
to roll up following various schedules (Ref. 5) 

b) A model in which the wake is represented by spiral 
vortex filaments (Ref. 6) 

c) A simplified model in which the wake is represented 
by cither vortex rings or doubly Infinite line 
vortices (Ref. 7) 

Section III examines the roll up process of the tip vortex in greater 
detail and presents an alternate treatment of the entire near wake and its 
roll up process with suggestions for a more elaborate approach to the problem. 

Section IV briefly discusses far wake effects and shows the dependence 
of wake geometry on the persistence of circulation in the far wake. 

Section V presencs the background theory used In the development of 
the various wake models. 

Section VI summarizes conclusions from this study and presents 
recommendations for further investigations. 
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SECTION II VORTEX WAKE MODELS 

The various wake models and their evolution to the final version using a 
combination of vortex sheets and vortex filaments are discussed In detail In 
Rafs. 5 and 6 which form Volume II and III of the present aeries. The first 
wake model used was originally developed In Ref. 8 for the treatment of wind 
turbine wakes. The same model was applied to a hovering helicopter rotor 
and preliminary results using this program were presented In Refs. 9 and 10. 
Ref. 9 reported progress at midpoint In the original contract schedule, 
comparing results using the original program with the experimental data of 
Refs. 2, 11 and 12. Discrepancies were noted and plans for future research 
outlined. Ref. 10 reported progress at a later point in the contract 
schedule, and Included results from the first two of five free vortex 
models developed during the course of the Investigation. The vortex sheet 
models of Ref. 10 contained- an Inadvertent doubling of the' strength of the 
far wake which resulted in a small error, corrected in the later models dis- 
cussed in Ref. 5. Again, discrepancies with the experimental results were 
noted which were attributed to the termination of the vortex sheets In less 
than one spiral. Attempts to extend the solution beyond 360* were un- 
successful because of lack of convergence. The solution was based on 
Iteration of the wake geomet / which required an appreciable amount of CPU 
time due to the complexity of the distorting vortex sheet model. The problem 
with convergence was believed due to the roll up of the vortex sheet which 
was not being properly modelled because of the coarseness of the mesh re- 
quired to keep CPU time to manageable proportions, of the order of an hour. 

The four vortex sheet models developed in Volume II nay bo summarized 
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as follows: 

model 1 Essentially the sane as Ref. 8, In which the near and 
intermediate wakes are represented by a segmented vortex sheet 
connecting a tip and root vortex. The far wake is represented by 
vortex cylinders as in Ref. 7, 

MODEL 2 The near wake is treated as in Model 1 until an azimuth 
of ■ 255*, after which it is rolled up into three vortex filaments 
as in Ref. 7. The location of these filaments at ■ 255* is 
determined as if the filaments had rolled up immediately at ” 0. 

MODEL 3 Same as Model 2, except that the near wake is rolled up 
before first encounter, at <$> •■ 45*, and is assumed to remain in the 
rotor plane until roll up. The geometry of the rolled up vortices 
at -p ■ 45® is determined as in Model 2. This avoids the lengthy 
computations required to establish the geometry of the vortex sheets 
and hence is computationally much more efficient than Model 2. 

MODEL 4 Same as Model 3, except that the full vortex sheet model. 
Including tip and root vortices, is used to establish the geometry of 
the vortex filaments at roll up. 

In addition a filamentary wake model. Model 5, was developed as dis** 
cussed in Volume III, Ref. 6. 

MODEL 5 The filamentary near wake is allowed to roll up as in Model 2, 
but to - 70*, after which the centroid of the partially rolled up wake 
is used to describe the wake geometry. The roll up process is similar 
to that of the Model 2 vortex sheets shown in Fig. 4. 
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The continuing InvestigaClons discussed In Section III Indicate that 
the roll up with a more finely modelled near wake with trailing filaments 
located approximately every ,3Z span at the tip would be much faster than 
with the relatively coarse mesh size with a minimum spacing of 2Z of the span 
used with the vortex sheet model. The roll up predicted by the vortex sheet 
model at the first Intersection (180°) Is shown In Fig. 4. Although the 
start of roll up Is clearly evident. It is not nearly as tight as predicted 
by the finer mesh model of Section III as shown In Fig. 10. The location 
of the tip vortex close to but above the blade most probably caused the 
oscillatory divergence of the solution. 

In the meantime, a separate investigation had' been conducted In which 
attempts were made to develop a simplified model of the wake which could be 
used to aid in clarifying the physics of the problem. This Investigation 
resulted In the development of the analytical technique described in detail 
In Ref. 7, with later results appearing in Refs. 13 and lA. The simplified 
program is capable of predicting wake geometry and air loads In a few 
seconds of CPU time and with reasonable accuracy. Velocities are 
computed in the wake immediately behind each blade only, and these 
velocities are averaged in order to obtain wake displacements. The 
averaging process in the near wake was based on experience with the 
more complete solution (see page 8 and Fig. 7) and avoids the need for 
computation of the wake displacements during the roll up process. 

Some of the results obtained with the simplified model are shown 
in Fig. 2. The computed results are in reasonable agreement with test 
data as regards bound circulation distribution and geometry. Some dis- 
crepancies exist in the region of the drop off of circulation near the 


ti'-aatifa 
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80% span point for reasons which are not clear but which are believed to 
be associated with the spanwise flow and real fluids effects discussed in 
Ref. 14 and in Section III. 




The simple model was based on the assiimption that the trailing wake \ U 


of the blade rolled up almost immediately according to the Betz criteria 
of conservation of angular and linear momenta (Ref. 7) and only the rolle 
up wake need be considered in the intenacdlate and far wakes. The near 
wake was treated as a relatively fine series of semi-infinite vortex fila- 
ments extending from the blade and was used only in determining the bound 
circulation distribution on the blade Itself in the presence of Induced 
velocities from the rolled up vortices. Fig. 5 Indicates schematically 
the simplified wake model. The velocities Induced by the vortex rings and 
cylinders could be expressed as elliptical integrals amenable to rapid 
series solution, consequently CPU time was kept to a minimum. 

In view of the success of the simplified model it was decided to 
use a similar roll up schedule for the more detailed models. The vortex 
sheets were therefore rolled up, again according to the Betz criteria, at 
various points in the near wake thereby avoiding the problems with the slow 
roll up predicted by the coarse meshed model discussed above. Volume II 
(Ref. 5) contains a detailed discussion of the various attempts to obtain 
satisfactory results with this model. Typical final results (Fig.. 6) 
show reasonable agreement with test data. 

It became evident in the development of this revised vortex sheet 
model that a problem existed in determining the local velocities con- 
tributing to the displacement of the tip vortex, either rolled up or in 
sheet form, between its generation by one blade and encounter with the 
following blade. The rolled up tip vortex is located at approximately 
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97. 5Z of the span and the computed velocities at the blade tip and at 95Z 
span typically are as shot^n in Fig. 7. One possibility would be to average 
the velocities at these two stations giving the results shown by the heavy 
line in Fig. 7. This was the basis for adopting the averaging process in 
the simplified solution of Ref. 7 and in Ref. 5 as a reascnable approximation 
£or determining the position of the rolled up tip vortex at first encounter. 
However such treatment of the tip vortex is an approximation only and not 
necessarily the true mechanism determining Its migration in the near wake. 

In an attempt to obtain a clearer understanding of these problems, 

Model 5 was developed as described In Ref. 6 with the wake represented by 
vortex filaments rather than vortex sheets. The original vortex sheet model 
was developed on the basis of experience with an earlier filamentary wake 
model developed in Ref. 15 for the forward flight case where it was found that 
a better representation of blade loads could be obtained If the wake at first 
encounter were represented by a vortex sheet rather than a concentrated curved 
vortex. However, in hovering flight, and at least for the two tladed case, 
the vortex at first encounter is sufficiently far below the blade so that the 
filamentary representation of the wake is adequate. Certainly the computational 
difficulties encountered in matching the edges of the segments of the dis- 
torting sheets are avoided. 

Results from the filamentary model are discussed in detail in Vol. Ill 
(Ref. 6), together with conclusions as to its limitations and that of the 
other models. The predicted wake geometry agrees well with the experimental 
data although the previously noted dip in bound circulation near 80Z span 
persists and requires somewhat arbitrary assumptions as to core size to 
Improve agreement with test data. 

It should be noted that the change in lift at the 80% span may be 
sufficiently rapid to trigger separation of the type first noted in Ref. 16, 


; 1 L ' i i;i h'-tf t r sr: ■. . ■ ■ r. >- v - i," 
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and discussed in Ref. U. when a close blodc/vortcx encounter occurs. 
Furthensore. there exists at the 80: span a wake induced radial velocity of 
the order of 105 of the tangential velocity due to blade rotation. The flow 
will therefore approach the blade at an angle of the order of 5* wh.'ch will 
introduce further second order effects siallar to the effects of sweep. A 
more coopleto lifting surface solution of the type discussed in Ref. 17. 
rather than a lifting line solution, is required in order to investigate 


such effects. 
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SECTION III VORTEX HOLL-UP 

. i 

Recent Interest in the fomatlon of \*ortlcos from fixed wind aircraft 

has stlnulated the developaent of coaputatlon techniques for predicting the 

geoQoCrles of rolled up three dlaonslonal vortex sheets. In these techniques 

the sheet is treated as a collection of line or point vortices and their 

notions tracked as they distort under their autual Interferenco velocities 

18 

using various coaputational nethuds. Sone of the earliest work is quoted 
in Ref. 19, pp. SS9-590. Subsequent efforts to duplicate these results with 
a finer grid resulted in chaotic notions, particularly in the tighter portions 
of the spiral. Introduction of artificial viscosity resulted In noro ordered 
solutions but, unfortunately, dependent on the degree of viscosity Introduced. 

Most recent efforts have concentrated on the EuIer**Lagrange solutions 
in which a series of point vortices are tracked in a Lagrangian frane of 
reference and referred back to on Eulerian fraoe for solution of the equations 
of flow. This **cloud In coll" technique was applied to the roll up of 
vortex sheets froa a wing In Ref. 20 and core recently In Ref. 21. In Ref. 22 
the two dlnenslonal approximation to the rotor wake developed in Ref. 7 was 
applied to the computation of the roll up of an assumed straight wake from a 
rotor blade as it descended in hovering flight. A typical solution Is shown 
In Fig. S which clearly indicates the roll up process. Time did not permit 
completing Che study to a converged solution. This work lx being continued 
using Che throe dimensional model of Ref. 7 which allows for curvature in the 
wake. Preliminary work has concentrated on determining the effects of coll 
size as one measure of the degree of artificial viscosity introduced in the 
solution. Dependence of the solution on the assumed cell size remains a matter 


of concern. 
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In chis section we will discuss a further approach to the probica of 
vortex sheet roll up In which the wake is i&odollcd as a curved scries of 
vortex flXaaents and their roll up predicted using the Biot-Savart relation- 
ships. In Ref. 6 the conplcte solution has been obtained for a series of 
spiral vortex filasents, but only a liaited nuaber of such filoacnts could 
be used in order to keep the CPU cine reasonable. In this section the roll up 

of the near wake only will be considered and a sisplifled approacit taken in order 

to nlnlnise the required conputer time and clarify the nature or the problem. 

Referring to Fig. V, ic is necessary to compute interferer '.e velocities 
between any two vortex filaments everywhere in the first spiral and integrate 
these velocities to obtain displaceaenc. If Ic is desired to compute the 

velocity Induced at any point A by another vortex filament B, tnen it is is 

clear that the velocity at A due to B would be primarily due to that portion 
of B closest to A. The rest of the spiral vortex filament B nay be approxi- 
mated by a vortex ring and the velocities computed readily by the teciinlques 
discussed in Section V using the Blot-Savart relationships and logarithmic 
series solutions for the resulting elliptic Integrals. The displacement of 
any vortex nay be computed by integrating the total induced velocity on a 
vortex due to contributions from all ocher vortices in the near wake over an 
increment of tine represented by a small change In azimuth, starting from the 
blade in question, using standard techniques of integration such as fourth 
order Runge/Kucta formulas. With this simplified tsodcl it is possible to set 
up a much finer wake structure, thereby presumably achieving a more rculisClc 
wake roll up. 

This technique was first used to examine the roll up of the curved 
vortex generated from the tip of the blade using a fine distribution of 
filaments. A typical result is shovm in Fig. 10. Nineteen filaments were 
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generated froa the outer alx percent of the blade with an assuaed core size 
of IZ of blade span. Most of the vortlclty Is contained in the nuabered 
vortex filaments and in particular the first six as is evident from the table 
showing the strength of each vortex fllaaent, g, where g is equal to V/QK^ 
with a aoxlaua value of .02 at 94X span. It cay be deduced froa Fig, 10 
that the vortex rolls up very rapidly, first rising and then descending. 

Roll up apparently occurs a few chord lengths behind the rotor blade, as cay 
be expected froa the experiaeatal evidence. 

The vertical dlsplaceacnt of what is apparently the vortex core is 
characteristic of a curved vortex sheet and not a straight sheet. In the latter 
case, it cay be shown that no vertical displaccacnt of the centroid of vor- 
ticity will occur if the effect of blade bound vorticity on the wake is 
neglected. 

The wake displacecents due to the bound vorticity of one blade will, 
to first order, be cancelled by that of the following blade, resulting in no 
net dlsplacecent of the vortex fllaaents, whether curved or straight, due to 
bound circulation. Fig. 11 shows the dlsplacecent of a free vortex with 
.and without the presence of bound vorticity coeputed using the analytical 
techniques discussed in Section V. The net effects are evidently scall 
everywhere in the wake and essentially zero at the first encounter with the 
following blade. Consequently the bound circulation has not been included in 
the results of Fig. 10 in order to clarify the core Inportant effects due to 
the mutually induced velocities of the free vorticities. It should be noted 
that Che bound circulation is Included in all Che coaplece solutions dis- 
cussed in Vols. II and III. 

Two factors renain of concern in the roll up predictions: the step size 
used in the integration of velocity in order to obtain dlsplacecent and the 
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aasuned vortex core size. The effect of core size on the dlsplaccaent of the 
centroid of vortlclty ia shoun in Fig. 12. Evidently the diaplecenent ia 
dependent on the aasuned core aize. When zero core size ia used the results 
tend to becone chaotic near the centroid of the spiral as previous investi- 
gators have found. Furthensore the solution beconea nore sensitive to the 
interval size used in integration whereas this ia not the case for a solution 
containing a sufficient amount of artificial viscosity as shown in Fig. 12. 

For a two bladed rotor Fig. 2 indicates that the tip vortex at first 
encounter is located approximately 5 to 6Z of the span below the following 
blade, as shown both analytically (Ref. 7) and experimentally (Ref. 4). 
Consequently, the differences in the vertical displacement of the centroid 
of the rolled up vortex with various assumptions as to core size, of the 
order of less chan IS of the blade span at first encounter, taay not appear 
to be crucial. However blade airloads arc sensitive Co the location of the 
tip vortex at this first encounter and the dependency of the solution on 
assumption as to core size remains of some concern. Experimental evidence 
(Ref. 23) indicates vortex core sizes of the order of 12 of the span, but 
these measurements are presumably of the rolled up vortex core, which nay 
not necessarily be Che core size of the filaments which model the sheet 
as it leaves the blade. A great deal more analytical and experimental 
investigation is necessary before this problem nay be completely resolved. 

It is possible that simplified forms of Che Kavier-Stcites equations for 
spiraling vortices will have to be developed in order Co obtain a better 
understanding of the phenomenon of roll up and migration of the tip 
vortex generated by a rotating blade. 

As a first attempt to extend this investigation Co the complete wake, 
an approximate solution was obtained using 24 vortex filaments for the near 
wake. A converged solution was first obtained using the fast free wake 
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tuchnlques of Ref. 7, whoro tho naac wake is asausied to coll up alooat 
instantly. The resulting wake geoaetry was then used to predict tho 
intecBedlato and far wake velocities at the blade and these were added 
to those computed froa the fine (now curved) near uaka by interpolation. 
Tho results are shown In Fig. 13 which clearly shows tho roll up of tho 
tip vortex and the Initial roll up of the first center vortex. Tho 
sudden Jump in displacement between the tip vortex and the center sheet 
has frequently been observed in experlnental Investigation. 

The next step will Involve modifying either the fast free wake modal 
of Ref. 7 or tho complete model of Refs. 5 or 6 to include tho fine near 
wake characteristics indicated in Fig. 12 and completing tlie iteration for 


blade loads. 
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SECTXO!; IV FAR W^VRE EFFECTS 

Another effect requiring ooro dotallod modelling than la possible 
with the programs of Refs. 5 and 6 is that duo to tha far wake. In the 
free wake analyaos, the far wake io represented by acmi^inf inite vortex 
cylinders whoso contributions to the Induced velocities at the blade are 
small. However tha influence of these vortox cylinders io more noticeable 
on the geometry of the intermodlato wake and may thoreby indirectly affect 
the induced velocities at the blade. Also it has been pointed out in Ref. 3 
that experimental ovidence indicates a posoib.le expansion of the wake 
several spirals below the rotor. 

The fast free wake technique of Ref. 7 allows a closer examination of 
the intormadiate wake by the introduction of many more spirals than is 
foasible with the conplote solutions. Fig. 14 shows the standard solution 
using the vortex cylinder representation for the far wake for the usual two 
spirals per blade before starting the for wake and for a case with nino 
spirals which may be expected to bring the far wake sufficiently for away 
from tho blade so as to make its effects negligible. The results chow that 
the solution is esoontially independent of the number of spirals, providing 
at least four arc taken in the intermediate wake. Fig. 15 shows the some 
solution but without tha constraint of a fixed far wake. The effect on rotor 
performance is small as may be expected, but it is interesting to note an 
expansion of the intermediate wake which may bo attributed to the Intormingling 
of the outboard center and tip vortices In a typical Kelvin /Ho Imho Its Insta- 
bility. Continuation of tho solution with core spirals may bo expected to 
show an oscillatory expansion and contraction. Further experimental and 
analytical investigation of the growth and dissipation of the far wake would 
appear to bo desirable. 
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OF P0C3 QUALrrv 


SECTIOM V DEVELOPMENT OF AKALYTICAL TREATMENT 


Velocity Induced by a Vortex Ring 

ha discussed In Section III the roll up of the near wake oay be com** 
puted by considering the velocities on one element, A, (Fig. 9) of the near 
wake due to all other elements by treating the spiral wake as a series of 
vortex rings. This assumption neglects the effects of distortions In those 
portions of the wake far from the element In question, evidently a reasonable 
approximation since these distortions are small relative to the radius of 
the spiral. The relative positions of all elements in the wake in the 
vicinity of the element A, however, arc correctly modelled. 

If r Is the radius of the vortex ring B of Fig. 9, n Is the radial 
location of A and z Is the relative vertical position assumed Independent of 
(no distortion of the ring) then the vertical component of velocity Induced 
at n due to the vortex ring B Is (see, for example. Ref. 24) 



fZit 

r(r " n cos ») 

Q (n^ + + 2 ^ - 2rn 



With the substitution ■it - 2 'J; such that COS 0*2 sin^ - 

“”Tr/2 « 

T - -r(r + n - 2n s1n -u) 6n> 

h ^ — 9 — b ^ 

Jtt^ 2 n +2 - 2rn - 4rn sin ii, 
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ORIGINAL PAGE IS 
‘ OF POOR QUALfTY 

and since the Integrand is ayraactrical about 0 « 0 the Integral need be 
evaluated only from 0 to tt/2. After some manipulation one obtaina 
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These integrals may be evaluated by standard methods for complete 
elliptic intervals of the first kind 
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which, after substitution for A and B becomes 



rn 


t K - E[1 - .5k^0 + r/T^)l/(l - k^)l 


from which w may be evaluated for any value of k. 

In the absence of a computer with a library of elliptic functions it 

is convenient to use scries solutions for E and K. The cost suitable for 

the present application are Cayley’s logarithmic series because of their 

2 

rapid convergence, although care must be taken in the vicinity of k — >1. 
From Ref. 25, 777.3 and 777. A, and defining 


F ** In (- 


I - k‘ 


), 




E . 1 + .5(F - .5)(1 - k^) + ^ (F - 1 - \^)0 - k ) 


K . F + .25(F - 1)0 - k^) . (F - 1 - 1)0 - k^)^ + 


The radial component of Induced velocity at n is 




z r COS (|> d6 


(n^ + - 2m cos ) 


3/2 


which, following the same procedure used above may be evaluated as 


“ " 4SF • 2^/^ tE(2 - k^)/0 - k-=) - 2C, 
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A self induced velocity will exist on the -curved vortex, which 

A 

Ref. 26, may be expressed as 

" 7^ C^n(8r/Ar) - J-] 

where Ar is the core radius. Ref. 27 suggests a value of ~ of .025 
as typical for helicopter rotor loadings. The experimental results of Ref. 21 
indicate that could vary from .01 to .03, depending on the proximity of 
the vortex to first encounter, hence the term in brackets could vary from 
about 5.5 to 6.2. Since the singularity is logarithmic, and the contribution 
of the self induced velocities is small, the solution is not sensitive to an 
exact definition of ~ . A value of WS ■ [6] has been used for all 

the cases of Fig. 12. 

In Ref. 7 it was shown that the drop off of bound circulation at the 
tip of a rotor blade nay be closely approximated by an expression of the 
form 

. 0 . , 2 ) 1/3 

0 

which has been used in determining the strength of the trailing vortex 
system for the computation of the tip vortex roll up of Fig. 10. However 
in the case of Fig. 13, the computed bound circulation distribution has been 


used 
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Effect of Bound Circulation 

In order to examine possible effects of the bound circulation on the 

near wake displacements, a simple case of two parallel infinite blades of 

chord b operating in a wind stream of velocity V and subjected to a uniform 

downwash u will be considered, but with the blades treated as airfoils, or 
0 

lifting surfaces, rather than lifting lines, since the effects of the bound 
circulation are expected to be of importance primarily in the vicinity of the 
surfaces themselves. 

The bound circulation may be represented by a series (see, for example. 
Ref. 19,), of the form 

ft ® 

y(x) » A cot I- + y sin n9 
' 0 2 n 

^ n“l 

where X = ^(1 - COS 8). 

The velocity induced on the airfoil by the element of vorticity y(x)cIx Is 

A A 

v(x) = J Y cos ne 

For a flat airfoil in steady flow, only A is nonzero. 

0 

It is desired to find the dlsplacemeni: of any element in the wake 
between the airfoils when both are at an angle of attack a. 

Boundary conditions on the blade require that , for o. 


v(x) = Va 
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everywhere on the blade whence 

A = 2Va 
0 

. - sin 9 

and, since COu 0 ® -j . 0 

/ A o\t sin e 
y(x) - 2Va 1 _ cos 9 

a 2Va/ b - X 

"nr* 

Referring to Fig. 16, the velocity at any distance K behind the first airfoil, 
and p - C ahead of the second airfoil, due to an elenent of bound vorticity 

Y(x) on each airfoil is 

y(x)dx r g b - X P ~ ^ ^ 

^ [ (5 + b - x)" + (P - 5 + + z^ 

If a steady uniform velocity u^ exists perpendicular to the airfoils, for 
example from the rest of the wake, then the vertical velocity at 5 is 





. ^ - V 

or, since dt ^ 


it 


whence 

f x-b ^ 

^ ” I II 

5«0( x«0 


b - X 


b + g - X 


L (b + ? - x)" + z 


P - 

^ (P - C 


+ x)^ + z^J 
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which may be readily programmed fcr direct numerical integration. Fig. 11 
shows the results obtained for the case of a .1 radians on both airfoils. 
Evidently the displacements are everywhere small, and essentially zero at 
the second airfoil. Fig. 11c shows the displacements when u/V = has 
a typical value of .02. 

The actual conditions could of course differ from those of the 
simplified model used here. The bound circulation varies along the blade 
and particularly near the tip. The trailed tip vortex moves rapidly 
inboard and therefore may encounter a slightly different velocity field due 
to the bound circulation near the second blade. In order to examine the 
extreme of such a condition Fig. 11c shows the case where the bound circu>^ 
latlon on the second airfoil was set equal to zero. The solution then 
reverts to the familiar case of the downvash behind a llfing airfoil. The 
z displacements are still small, of the order of 20Z of the vertical dis- 
placement at first encounter fcr a two bladed rotor. The bound circulation 
has been included in the models used in Volumes II and III, although its 
effect was found to be negligible. Its neglect in the simplified solutions 
of Section III appears to be a reasonable approximation. 
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SECTION V! 

Cone Iu9 tons 

1) Free wake techniques of varying coaplexlty have been developed 
for predicting airloads on rotor blades* all of which agree well with 
experincntal observation of both the bound circulation distribution and 
wake geoaetry. 

2) Simplified Tcodels in which velocities are conputed only in the wake 
directly behind each blade and averaged to obtain wake displacements give 
results In agreenent with the core conpletc solutions in which wako velocities 
are computed everywhere In the near and intonsedlate wakes. 

3) Techniques liave been developed for predicting the migration of the 
tip vortex in the near wake for the complete solution, however uncertainties 
as to the exact roll up oechanisn and near wake geometry still remain. 

4) An analysis of the roll up acchanisn in the near wake using a 
simplified model which permits use of a fine trailing wake indicates that 
roll up of the tip vortex is essentially completed in a few chord lengths 
behind a blade. 

5) The vertical displacement of the centroid of vortlclty in the near 
wake depends on the assumed core size (artificial viscosity) of thr; vortex 
filaments representing this wake and the solutions became chaotic for small 
or zero values of core size. 

6) Use of a rigid far wake has little effect on the predicted blade 
loads but does affect the geometry of the intermediate wake. 

7) Effects of bound vortlclty on near wake diaplacemcnta are neg^lfilhlc 
and its neglect in the simpLifled solution is Justified. 
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Recogacndaciong 

1) Extend either the simplified or more detailed solutions to Include 
a sore rcaploto representation of the near wake vortex roll up and Itarata 
for a final converged solution using this taore complete wake geometry. 

2) Develop cx>re fortial techniques for including the effects of viscosity 
In the deteralnatlon of wake goonetry. 

3) Extend the solution to include a lifting surface representation of 
the blade for use with rotors having four and tware blades and In order to 
allow for the effects of apanwlso flow along the blade. 

4) Obtain cxperlaental verification of blade bound vortlclty diatrl** 
butlon and wake geometry for rotors with four and oore blades. 
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Fig. 2 Blade bound vortex distribution and wake geometry as predicted 
by methods of Ref. 7 compared with experimental results of 
Refs. 4 and 11. 
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Fig. 5 Geometry of model using vortex rings and cylinders to represent 
the wake. 

a) Side view of rotor wake model showing intermediate and far 
wakes formed from vortex spiral - 2 blades. Tip vortex only 
shown. 

■ Blade One Blade Two 

b) Plan view showing near wake 

c) Formation of intermediate wake 






cAfcJfcC 













ORIGINAt PAUL 
OF POOR QUALITY 


37 


t • • *.tn 

» 4. lima / f-i 


• M IM» 

aoLl. Of fofoo 'i» • »af«s 

M K.lSr •'» *•** 

f* 4# 

I t 4 7.J44* C-> 

t • ».3ii«a i-i 
] a.ta) • i.Tiaca s-l 
4 a.M • i.aiTri ao 
s $.tOf o 1.1333* C-1 
a a.iea • f.isrtr s>« 

• a.*a • s.iaMa (->4 
a a. *73 « t.iiiJt (•* 

• *.*71 • 3.:o*37 t-4 
a a. *7 4 *.M7«7 1*4 

• *.*47 4 «.«t*34 ■-* 

« a.tia • ].irs7 c*4 

• *.*1 4 a.«*77a a-4 
■ *.*37 4 4.1*13 1*4 
» *.*J3 4 1.7751* l-» 

4 a. *3 4 I.I3IS3 (-4 
« *.*47 4 1.44ft** 1-* 

• *.*43 « 4.41)43 k-9 
4 *.*4 4 4.1*34* 1-1 


f - 0 

'1* l***^«*»<****»^^“^3*< Jl- 1 


yC 



/• 




~} f - la* 

t 

1 

**>• .*•' 





/ 








A 



1 11 * i.:t 


■ a.*ti:aa y r. 

4)471 1-3 


/" ^ 


/ 


( r 


a < /7« 


1 \ Oi 


' 


\ 




> 


- A5* 


« 90* 


\ 


Fig. 10 Predicted roll up of tip vortex - near vake only considered 
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Fig. 11 Effect of bound circulation on vortex aigration between two 
lifting surfaces. 

(a) Effects of bound circulation only fron both surfaces 

(b) Sane but with bound circulation on originating surface only 

(c) Sane as (a) but including a steady inflow of X “ .02 

(d) Effect of steady inflow X ■ .02 without bound circulation 
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Fig. 12 Effect of core sice and interval size on migration of tip vortex 
centroid in near wake, $ « 0 to 180“. 

(a) Effect of core size 

X .02ft o .OlR 0 . 002:1 

(b) Effect of inter/al sire 

o of .05 tt Core size of .OlR 

AV of .025“ Core size of .OlR 

□ AY of ,05^ Core size of .002R 

7 AY of .025“ Core size of .002R 

(c) Effect of interval sire for core size of zero 

□ AY - .05t 7 AY - .025^ 
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Fig. 14 Effect of extent of intertaedlate wake with far wake 
2 spirals « .00456 9 spirals “ .00454 
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Fig. IS Effect of ellslnatlng far wake with extended Interaedlatc wake 
9 gpiralo ■ .0044 





